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SUMMARY 


Simple formulas are given for the stagnation-point convective heat 
loads in lunar return for two operational modes. The two modes of 
operation analyzed are typical of moderate heating and of nearly minimum 
heat loads, respectively. The values of the parameters in a simple two- 
parameter formula for the total-heat load are given in the lift-drag- 
ratio range of 0.2 to 1.0 and in the peak loading range of 2 g to lOg. 

For vehicles having a lift-drag ratio near 0 . 5 > which is typical of many 
proposed lunar return vehicles, the nominal mode had about 20 percent 
more absorption than the nearly minimum mode. 


INTRODUCTION 


In the course of preliminary trajectory analysis of the lunar 
return problem it was found possible to compress a large amount of 
stagnation-point heating data into simple formulas for a range of 
vehicle lift-drag ratios and a range of peak loadings which are of 
current interest for manned vehicles. Several modes of operation which 
span a range of heat loads are considered. 


SYMBOLS 


A reference area, sq ft 

c satellite speed at altitude of 50 miles, 25,800 ft/sec 

Cjj drag coefficient 

C^ constant in heating-rate formula 


D 


drag 



2 



g acceleration due to gravity, 32.2 :ft/sec^ 

h altitude, ft 

L lift 

n exponent in heating- rate formula 

Q heat absorbed per unit area at sta ( piation point, Btu/sq ft 

q heating rate per unit area at stagnation point, Btu/sq ft-sec 

R resultant force 

r radius of curvature, ft 

s longitudinal range, ft 

t time, sec 



V airspeed, ft/sec 

W weight, lb 

3" 1 scale height, 24,500 ft 

7 flight-path angle, deg 

6 - cJS? lb / cu ft 

p air density, slugs/cu ft 

0 roll angle, deg 
Subscripts : 

1 start of constant- altitude segment 


L 

1 

6 

1 

5 


1 


start of constant-G segment 
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2 end of constant-altitude segment 

2' end of constant-G segment 

Dot over symbol indicates differentiation with respect to time. 

PARAMETER RANGES 


The parameter range considered is indicated in figure 1 in terms 
of Chapman's corridor width. (See ref. 1.) Corridor width is a meas- 
ure of the error in perigee height of the approach orbit which is per- 
missible for a given load limit. Fixed angle of attack and fixed 
geometry are assumed in figure 1. Consideration of guidance accuracy, 
atmospheric density variations, and piloting errors will fix the cor- 
ridor width. Since the gross weight of the vehicle tends to increase 
with (L/D) max capability, the lightest vehicle will have the lowest 

(L/D) max capability consistent with the constraints of load limits, 
corridor width, and maneuver capability. If a load limit of lOg and 
a corridor width of 50 miles are adopted as extreme values, the high 
values of L/D need not be considered. Accordingly, attention has 
been restricted to values of L/D less than one. Simple formulas were 
not found below L/D =0.2 so smaller values than 0.2 are not included. 
If a lower load limit such as 5g is adopted (see fig. l), it should be 
noted that some minimum value of (L/D) max is required to meet the load 

limit even with perfect guidance (zero corridor). 

In short, small corridor widths allow low values of L/D, and 
small load limits force high values of L/D. Corridor width is inde- 
pendent of W/Cj^A so vehicle weight does not appear in figure 1. The 
parameter range of 25 = w/Cj)A 1? 400 lb/sq ft was adopted in the cal- 
culations to span a wide range of weights. 


MODES OF OPERATION 


The aerodynamic coefficients are not allowed to vary in the modes 
of operation considered. The values of L/D and 5 are thus fixed 
through each entry. Control is exercised by rolling the vehicle (and 
thus the lift) about the velocity vector. In general, a discontinuous 
jump in roll angle followed by a continuous variation is required. 
Formulas defining the required roll-angle variations are given in 
appendixes A and B. 
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Constant G 

A convenient lower bound for the total-heat load to be experienced 
in lunar return is provided by what will be called a constant-G 
(dG = 0) entry. In this entry (see fig. 2) the vehicle makes a constant 
roll-angle (d0 = 0) pull-up to the point of peak G. At this point a 
sudden roll is presumed which is followed by a continuous roll in such a 
manner as to keep the aerodynamic loading constant. The constant-G entry 
maintains during the slow-up the highest Reynolds number consistent with 


the load limit and hence the lowest ratio of friction to pressure drag. 

If the constant-G flight is presumed to continue to zero velocity, a L 
lower bound to the heat absorption is established. The speed to which 1 
constant G can actually be maintained (see appendix A) is near enough 6 
to zero in the parameter range of interest that the bound obtained can- 1 
not differ greatly from the sharper bound obtained in an entry with a 5 
dynamically possible low-speed portion. 


Constant h 

A second mode considered is the constant-altitude (dh = 0) slow-up. 
This mode is believed to yield heat loads more typical of actual one- 
pass entries than the quasi- limiting case of constant loading. In the 
constant- altitude mode a constant roll-angle (a0 = o) pull-up is fol- 
lowed by the constant-altitude slow-up, which is followed by a smooth 
glide to Impact. (See appendix B. ) As in tie previous case, a con- 
tinuous roll is presumed as the control technique. The mode Is char- 
acterized by sustained but steadily decreasing loads through most of 
the entry. 


Constant 0 

A third mode considered was the fixed-control (d0 = 0) mode. 

Results are given only for L/D =0.5 which is typical of many proposed 
lunar return vehicles. This mode is objectionable for reasons other 
than higher heating, but it is introduced to provide a measure for a 
relatively high total-heating mode somewhat ns an upper bound to comple- 
ment the lower bound provided by the constant-G mode. Of course, the 
constant-0 mode is not an absolute upper bound. 


POWER LAWS 


Examination of the trajectory data obtained in the parameter range 
of figure 1 has indicated that, for the modes of operation considered, 
the peak heating rate at the stagnation point during entry is 
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W= K qV6G y (la) 

and the total heat absorbed at the stagnation point during entry is 

Q = K q ^&G" X (lb) 


for 5 = w/c j)Ar. The values of the constants Kq, K^, x, and y 

were obtained from plots similar to the sample shown as figure 3* The 
scatter of the data is typical of that observed in the range, 

0.2 ^ L/D ^ 1.0, and in the range, 25 = 8r ^ 4-00 lb/sq ft. Just as in 
figure 3 the slope of the line fixes x and the vertical position 
fixes Kq, so do the slope and vertical position fix the exponents and 

constants in all cases. The values of Kq, K^, x, and y depend, of 

course, on the mode of operation and the operating L/D. These values 
are shown in figures 4 and 5 for two modes of operation. In reference 1 
somewhat more complex relations than those of equations (la) and (lb) 
are suggested, and the corresponding constants are evaluated for a num- 
ber of modes of operation. 

Much of the data was taken from numerical integrations performed 
by an IBM 7090 electronic data processing system. Some analytical cal- 
culations were made as outlined in appendixes A and B. The stagnation- 
point, convective-heating- rate formula used throughout was of the type 
used in reference 2 and is as follows: 


q = 17, 300 




3.15 


Btu/sq ft-sec 


\10 H 


where p is in slugs/cu ft, V is in feet per second, and a radius 
of curvature of 1 foot is assumed. The assumed initial conditions in 
every case were V = 36,500 ft/sec at h = 400,000 ft in a static, 
radially symmetric ARDC 1959 atmosphere. (See ref. 3-) la the analyt- 
ical calculations, an exponential atmospheric density variation with 
altitude is used. 


RESULTS AND DISCUSSION OF TYPICAL VEHICLE 


Equations (la) and (lb) yield an interesting result when applied 
to vehicles having a lift-drag ratio of 0.5 which is typical of many 
proposed lunar return vehicles. Total-heat loads at the stagnation 
point of these vehicles are shown in figure 6 for a range of G^ax 
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from 2 to 10. Only entries for positive lifu in the pull-up are shown 
in figure 6. Limiting cases for negative lift in the pull-up are shown 
subsequently. 

With the assumption that the values of q/\/& for the constant- 
altitude mode are representative of likely levels of heat absorption, 
figure 6 shows a rather small improvement to be found by resorting to 
the constant-G mode and a similar small penalty for the change to the 
constant-0 mode. In each case the change in heat absorption is about 
20 percent. 

The range in heat loads for a given mode is of interest when 
entries high and low in the entry corridor are considered. Figure 6 
represents undershot entries in the sense that for each level of G max 

the steepest entry has been made within the constraint of constant coef- 
ficients. The steepest entry is assured because upward lift is used in 
the entire pull-up. The overshoot boundary nay be defined by an entry 
in which use of downward lift through the initial pull-up just insures 
that the vehicle will not skip. After this Marginal pull-up the vehicle 
may maneuver to lower altitudes and make sin upright pull-up followed by 
one of the previously discussed modes of operation. An illustrative 
case is shown in figure 7 again for the constant-h and constant-G modes 
at L/D = 0.55* Some description of the complicated modes of operation 
corresponding to the overshoot boundary curves in figure 7 is necessary. 
The vehicle enters the atmosphere with the lift vector in a vertical 
plane pointing downward. The entry is so shallow that when the vehicle 
levels out there is just enough lift to overcome the resultant of the 
centrifugal force and the weight. Since this resultant is decreasing 
faster than the lift, the vehicle bites deeper into the atmosphere in an 
inverted attitude. By rolling to a normal upright attitude at increas- 
ingly later times the vehicle may make pull-ups at increasingly higher 
values of G followed by the usual slow-up phases of the constant-h 
and constant-G modes. The heat absorption indicated as overshoot 
values in figure 7 corresponds to entries mac.e in this manner. It is 
apparent from the figure that the penalty for overshot entries within 
a given mode increases with the loading level of the entry, reaching 
appreciable values (larger than 20 percent) at the higher- loading 
levels. It was not found possible to accural .ely represent entries of 
the type considered in figure 7 by formulas similar to equations (la) 
and ( lb ) . 

Of course, it should be remembered that all results refer to the 
convective component of the heat load. Inclusion of the radiative 
component would add to the values shown. 


t-3 i — I VO < — I LTV 
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CONCLUDING REMARKS 


Simple relations have been given for the stagnation-point convective 
heat loads in lunar return for two operational modes • About 20 percent 
more heat absorption was found for a nominal mode typical of one-pass 
sustained loading entries (constant altitude) than for a mode with nearly 
minimum absorption (constant loading). For vehicles having a lift-drag 
ratio near 0.5> operation with fixed controls resulted in total heat 
loads 20 percent higher than those experienced in the constant-altitude 

L mode . 

1 
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1 Langley Research Center, 

5 National Aeronautics and Space Administration, 

Langley Field, Va. , June 9> 19&1. 
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APPENDIX A 
CONSTANT- G MODE 


The limiting case of constant G, so-called, can he easily handled 
analytically for the constant-trim, variable-roll mode of control. 

Constant trim implies nearly constant aerodynamic coefficients in the 

hypersonic speed range so that constant dynamic pressure is required 

for a constant- loading level. Thus, L 

1 

G ’ 5 = f * (d) Hf ‘ c,3n3tant < M > , 

The trajectory profile is sketched in figure 2. The condition that 
dG = 0 may be written as follows: 

^(pV 2 ) = pV 2 + 2pW - 0 (A2) 

Assuming an exponential atmosphere and neglecting gravity- induced 
speed changes gives 


P = -3pV sin 7 (A3a) 

-V = | g (A3b) 

Upon combining equation (A2) with (A 3 ), introducing s = V cos y, 
and integrating, the distance traveled on 12 ' (see fig. 2 ) is 
found as 


3s 


1 ' 2 ' 



1 + 



- 1 


(A4a) 


Equation (A4a) assumes the constant dec< leration to be maintained 
to zero velocity, which is incorrect but intj-oduces only a small error. 
For the large values of V,, which are of concern in lunar return, 


1 ' 2 ' 


Vl' 

2 §e 


(A4b) 
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which is just the result to be found by assuming 7=0 through the 
deceleration period. This result implies that 7 remains small through 
most of the arc 1 T 2*. 

The question of just how long the constant deceleration can actually 
be maintained is of interest. On 1 ' 2 ' the required flight-path angle 
can be found from equations (A2) and (A3) as 

(A 5 a) 


(A 5 b) 


Since 7 is small and remains small except in the neighborhood of 
point 2 1 } 


sln 7 - - jfc 0 


or as 


sin 7 = 


2gG u 


-2 


(§) £ 


cos 77 7 « - -§- p (A 6) 

08 r 

is the required rate of change of flight- path angle. Equation (A6) is 
assumed to be satisfied by rotation of the lift vector about the velocity 
vector in a slow continuous roll. The governing equation for 7 is 

— 7 = — cos 0 + (u 2 - l)cos 7 (A 7 ) 

8 W 

The value of 0 is to be varied in the manner required to satisfy 
equation (A 7 ) since it is the only quantity not fixed by the constraint 
of constant G. The condition cos 0 = -1 defines point. 2 ', the place 
beyond which constant G cannot be maintained. At point 2 * the lift 
vector is pointing downward in the vertical plane. Conditions at 
point 2 f are defined by combining equations (A3a), (A6), and (A 7 ) into 
the following relation: 

2 w tan 7 2 ' + W = ( U 2 * ■ 1 ) COS 7 2 ' ( A8a ) 

This transcendental equation for 7^, has supercircular and sub- 
circular solutions. The subcircular solution corresponds to point 2 1 . 
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A convenient approximate solution to equation (A 8 a) exists when 
and cos 7 ^, » 1 as 



“ — 

tan > , = i 

£ + -\/l + (lif 

2 ' 2 

D GV \D/ 


u§, « 1 
(A 8 b) 


Precise solutions to equation (A 8 a) are shown in figure 8 . The 
case of low values of L/D and high values of G will yield a measure 
of the smallest velocity range through which constant G can be main- 
tained. Figure 8 shows that, for L/D = 0.2 and Gj,^ = 10, - 7 2 , = 8.2° 

at u 2 , = 0.4. Even in this unfavorable case the fraction of satellite- 

speed kinetic energy remaining at point 2 ' is small. 


If stagnation-point total heat absorbed is calculated on the basis 
of uniform deceleration through the velocity range from u-, , to zero, 
then 1 


Q 


1 ' 2 ' 


O r v n 

-i\ 25_ V 

g \jD/W n 


(A9) 


when 

4 - 04 \/?v n 

The heat loads absorbed in the slow-up for the case of dG = 0 in 
figures 6 and 7 were computed with the use of equation (A 9 ) which ass ume s 
constant deceleration to zero velocity. The previous calculation of 
u 2 , shows that u|, is always small for L/D ^0.2 and G S 10. The 

conclusion to be drawn is that the values of Q in figure 6 are lower 
than those experienced in an entry with a dynamically possible low- 
speed portion. The values of Q yield a lower bound but not a greatest 
lower bound. However, the neglected contributions to Q appear to be 
small enough so that the lower bound provided by the constant-G mode 
is not far from a greatest lower bound for the constant-trim, rolling- 
control mode of operation. 


\J1 H ONP C 
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APPENDIX B 

CONSTANT-ALTITUDE MODE 


L 

1 

S' 

0 

1 

5 


The calculations of the characteristics of this mode of entry 
assume that the altitude control is a continuous roll about the velocity 
vector. The governing equation on the level flight portion is thus 


0 = — cos 0 + (u^ - l) (Bl) 

W 

which is the same as equation (A7) with y set to zero. A limiting 
condition in this mode of entry occurs when cos 0 = -1 is required at 
the bottom of the first pull-up. A necessary condition at the bottom 
of the first pull-up is the following inequality: 



(B2) 


For shallower entries than those for which (B2) is satisfied, a 
constant- altitude portion cannot be continued beyond point 1. In fact, 
there is a sharp but continuous transition in a narrow entry-angle range 
between high-apogee, long-range, skipping entries and those that satisfy 
inequality (B2). 

Although the constant-altitude mode is merely illustrative of one- 
pass sustained loading entries, the assumed rolling mode of control has 
a more general interest. No changes in longitudinal trim and hence no 
aerodynamic pitch controls are required. Of course, the changes in 
heading angle and the buildup of lateral displacement caused by the left- 
over horizontal component of the lift complicate the guidance problem. 


The elementary integrals of the motion along the constant-altitude 
arc 12 are as follows: 



(B3a) 


(B3b) 



12 


%2 


C q ^c n+1 uf 1 - u§+l 

(%« 


where 


(B3c) 



1 + 


L P 1 C 


D 25r 


L 

1 

6 

l 

5 


at the start of the smooth glide. 
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